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The phase transition kinetics in three phase systems was investigated using the numerically efficient cell dynamics method. A phase- 
field model with a simple analytical free energy and single order parameter was used to study the kinetics and the thermodynamics 
of a three-phase system. This free energy is able to achieve three phases coexistence, which for simplicity we call a, /3 and y phases. 
Our study focused on the kinetics of phase transition rather than the nucleation of a seed of a new phase that was introduced into 
the matrix of the old phase when the relative stability of the three phases were changed. We found dynamical as well as kinetically 
arrested static scenarios in the appearance of the macroscopic metastable phase. A few other interesting scenarios of the kinetics of 
phase transition in this three-phase system will be demonstrated and discussed. 
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■1. Introduction 

The study of the phase transition in a multi-phase system has 
been studied for many years. In particular, the multi-step phase 
transition which involves the formation of intermediate long- 
lived metastable third phase when the phase transition from one 
|phase to another occurs has been studied for more than a cen- 
tury |[li|2l- Recently, a renewed interest in the phase transi- 
tion in the multi-phase system lH, 0, S 0] has emerged. The 
formation of the thermodynamically metastable third phase is 
; important academically as well as industrially because many 
industrial products are in a long-lived metastable state |3]. 
■ In this report, we use the cell dynamics method employed 
previously to study the kinetics of the phase transition in a sim- 
ple three phase system |01 and study the kinetics for more wide 
spectrum of free energy landscape using phase-field model, 
which has been extensively used to study various scenarios of 
phase transition (sl |9l [lOl [llj, ll2[] . We employ cell dynamics 
method as we want to consider the evolution of multiple nucleus 
with general symmetry while the previous authors, instead, con- 
sidered the evolution of one-dimensional traveling wave solu- 
tion, m H [il [H, llTT of single nucleus. We pay special 
attention to the formation of metastable third phase during the 
phase transition between first to second phases. 

2. Model free energy for three-phase system 

In order to study the phase transition kinetics in the multi- 
phase system, we will use the partial differential equation called 
isothermal phase-field equation for the non-conserved order pa- 
rameter 111]: 

# _ 

'dt ~ ~~5^ 



(1) 



where ij/ is called phase-field and is actually the non-conserved 
order parameter and T is the free energy functional (grand po- 
tential), which is usually written as the square-gradient form: 



dr 



(2) 



where the local part of the free energy /(i^) will be specified 
later and will realize the multi-phase system. As a special solu- 
tion, Eq.([Tli is expected to allow an interface-controlled growth 
with nearly constant velocity l/19i] . It should be noted that these 
dynamics Eq. ([TJ always guarantee that the total free energy 
decreases monotonically because lEoll 

f^?.. = -fm^..O (3) 
dt J 6i(r dt J \Siffl 

by reducing the local free energy / or the surface tension that is 
proportional to (Vi^)^ in Eq. (|2|l. 

We will consider the simplest multi-phase system with only 
three phases and introduce a local free energy function f{i//) of 
triple-well form as a function of the non-conserved order pa- 
rameter ifr. We extend the analytical free energy /(ifr) proposed 
by Widom 121 (l in order to achieve the relative stability of three 
phases: 
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which can realizes the two-phase and three-phase system ac- 
cording to the magnitude of two parameters and Ay. We 
have three phases a with the free energy at (/tq, = - 1 , jS with 
at tf/p ~ 0, and y with fy at iffy = 1 . Now the free energy 
landscape consists of two wells a and y-wells, or three wells a, 
P and y-wells. The parameter A^ controls the relative stability 
of the intermediate fi phase while the parameter Ay controls the 
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stability of the y phase. Because the role of j and a is inter- 
changeable, we will only consider the cases when the y phase 
or the intermediate yS phase are the most stable thermodynamic 
phase. Therefore, we will restricted to Ay > as the free energy 
at the 7 phase is given by 



fy = /((A = 1) = -3 Ay 

while the free energy of the a phase is always fixed to /„ 
The free energy of yS phase is approximately given by 



1 2 

0) = -Ap - -Ay. 



(5) 
= 0. 

(6) 



Several typical shapes of the free energy function fiif/) for sev- 
eral sets of the parameters A/j and Ay are shown in Fig.[T] 
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Figure 1: The model triple-well free energy f{if/). In this figure, we show three 
cases: Case A (A;, = -0.3, Ay = 0.05)/„ < //j < fy. Case B {Ap = -0.2, 
Ay = 0.15) fa <fi3< fy, and Case C (hp = 0.2, Ay = 0.05) fa < fp < fy 

We observe from the figure [T] that there are typically three 
cases for the free energy landscape according to the relative 
stability of three states. Case A when fp < fy < fa is a rather 
special case as the intermediate [3 phase is most stable. Case 
B when fy < fp < fa has been considered by Bechhoefer et 
al. flS*] and Celestini and ten Bosch fTT] though they considered 
only a special traveling solution. Case C when fy < fa < fp i& 
the most interesting and relevant to thephase transition of sev- 
eral soft-condensed matter system f^, '?] and has been studied 
by Evans et al. 1.16. .17.1 for the kinetics of conserved order pa- 
rameter and by the present author ^ for the non-conserved or- 
der parameter when fa = fy The phase diagram of this system 
in the two parameter space (p, <5) is shown in Fig.|2] 

Similar triple-well potentials of different functional forms 
were used by several workers to study the nucleation 
and the metastable phase formation using phase-field model of 
non-conserve J? , 13 , 14 , Tsf] and conserved | TgI 12 ] order pa- 
rameter. 



3. Results of numerical simulation using cell dynamics 

3.1. Cell dynamics and phase-field model 



According to the cell dynamics method 112311 . the partial dif- 
ferential equation ([T]i is transformed into the finite difference 
equation in space and time: 
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Figure 2: The phase diagram for the model free energy Eq. gj. Case K:fp < 
fy < fa. Case B:/y < fp < fa, and Case C:/y < fa < fp- The triple point is at 
A^ = and Ay = 0, where all three phases coexist with fa = fp = fy = 0. 



where the time t is discrete integer and the space is also discrete 
and is expressed by the site index (integer) n. The mapping F 
is given by 

F[(A(f, «)] = gm, «)) + \ [« Ht, n) » -ipit, n)] (8) 

where the definition of << * >> for the two-dimensional square 
grid is given by 



(9) 



i=nn 



/=nnn 



with "nn" means the nearest neighbors and "nnn" the next- 
nearest neighbors of the square grid . Instead of the original map 
function g (i//) = i/r - 1.3 tanhi/r ll23ll . we used the map function 
that is directly derived from the free energy landscape /((/'): 



dtfr 



(10) 



+ l,n) = F[ip(t,n)] 



(7) 



which is essential in order to study the kinetics of phase transi- 
tion when a subtle balance of the relative stability of the three 
phases in the three-phase system plays a crucial role 10,0. Ell]. 

3.2. Kinetics of phase transition in a three-phase system 

Since we are most interested in the evolution or regression of 
the metastable phase during the phase transformation after nu- 
cleation, we will only examine the kinetics of phase transition 
when various composite nuclei made from three phases are pre- 
pared in the materials using the cell dynamics equation Eq. O 
and the model free energy defined by Eq. (|4]i (Fig.[TJ. 

3.2.1. CaseA(fp<fy < fa) 

We have incorporated the above free energy Eq. ^ into the 
cell-dynamics code Eq. O used previously to study the KJMA 
(Kolmogorov-Johnson-Mehl-Avrami) dynamics ll26ll . In this 
case, the intermediate middle y6-well of the free energy land- 
scape in Fig. [T]is most deep and the left and the right well is 
shallower In this case the /? phase is most stable and the y 
phase and the a phase are metastable. 

In figure|3] we prepare the metastable (white) y strip of width 
2 1 embedded in the middle of the metastable (black) a phase 
with size 128 x 128 and observe the evolution of the system. 
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Throughout this paper, we will use 128 x 128 system where a 
periodic boundary condition is imposed. We use the potential 
parameters - -0.2 and Ay = 0.15 which correspond to the 
curve for Case A in Fig.[T] Figure [3] shows that the most stable 
(gray) /3 phase appears spontaneously at the interface of two 
metastable y and a phases during the evolution. This stable /? 
slab continues to grow and invades the metastable y as well as 
a phase. Finally the whole material transforms into the stable /3 
phase. 




Figure 3: A Gray-level view of the time evolution of a three-layer stripes when 
= -0.2 and Ay = 0.15 (Case A). The black area is the metastable a phase, 
the gray area is the stable j8 phase, and the white area is the metastable y phase. 
Initially the less stable 7 phase is sandwiched by metastable a phase at f = 0. 
Even though, the initial state does not contain the stable /? phase, it appears 
spontaneously at the metastable a-y interface and continues to grow. 



3.2.2. Case B (fy < fp < fj 

Now the right y well is deepest, and the middle /3 well is the 
next and the left a well is most shallow (Fig.[T]|. This staircase 
configuration of the free energy landscape was previously stud- 
ied by Bechhoefer et al. 1 13] and Celestini and ten Bosch jlm . 

In Fig. |4]we start from the same stripe configuration as that 
in Fig. [3] Now the moving a-y interface unbind into a pair of 
moving a-/3 and moving y6-y interface and a macroscopic slab of 
metastable j8 phase appears. The small oscillation in the phase 
field 1/' in the y phase is probably due to the interference ef- 
fect because it is confined. Since the /3-y interfacial velocity 
is slower than the a-/3 velocity (Fig. |5]l, the metastable j8 slab 
appears and continues to grow. 



V 




Figure 4: The cross section view of the stripe structure similar to Fig. [3] The 
most stable y phase with i/< = 1 is sandwiched by the least stable a phase with 
i// = -I. In this case, dynamically stable /} phase with i/r =: appears even 
though it is thermodynamically metastable. 



Figure |5] shows the time evolution of the position of the a-/3 
interface defined hy iff - -0.5 and the /3-y interface defined by 
(A - 0.5. Two nearly straight lines with different slopes indicate 
that the two interfaces propagate with different constant veloc- 
ities. Since the /3-y free energy difference is smaller than the 



a-/i free energy difference (Fig.[T]), the y6-y interfacial velocity 
is slower than the a-yS interfacial velocity (Fig.|5]) as they are ex- 
pected to be proportional to the free energy diff'erence 1 191, 251. 
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Figure 5 : The time evolution of the position of a-/i interface and /J-y interface. 
The position is measured from the bottom of the 128 X 128 pixels of Fig. |5] 
Two nearly straight lines indicate the propagation with constant velocities that 
are given by the slopes. 



Our numerical simulation directly confirmed the existence of 
a dynamically stabilized metastable state predicted by the the- 
oretical calculations of Bechhoefer et al. 11311 and Celestini and 
ten Bosch |14]. However, the appearance of this dynamically 
stabilized metastable y6 phase is not due to the special symmetry 
of planer wave front assumed by those authors 111311 1411 . Figure 
|6]show the evolution of the multiple of circular nucleus of most 
stable y phase embedded in least stable a phase. Again, the 
metastable y6 layer appears around the growing y core sponta- 
neously and will grow. 




Figure 6: Evolution of the most stable circular nuclei of y phase (white) em- 
bedded in least stable a phase (black). Again the dynamically stabilized macro- 
scopic metastable /} layer (gray) appears around the y core and grows at the y-a 
interface. 



3.2.3. Case C(fy < fa < fp) 

In this case, the right y-well is deepest, while the left a-well 
is the next, and the middle /3 well is the shallowest and least 
stable (FiglTJ. This free energy configuration with metastable 
intermediate phase or buried metastable phase was previously 
studied by a few authors flSEB. 

Now the direct phase transition from the metastable right 
a phase to the stable left y phase (Fig. [T]) is prohibited from 
Eq. (O as the reaction path must go through the intermediate /3 
phase with higher free energy. As a consequence, a composite 
nucleus of a stripe made from the most stable y phase of width 
21 sandwiched by the next stable a phase of width 33 embed- 
ded in the center of least stable /3 phase as shown in Fig. |7] will 
be kinetically arrested as the most stable y phase cannot invade 
the a phase and cannot grow. Even though the a - y front is 
freeze, the a - /3 front can grow as the a phase can grow by 
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consuming the surrounding fi phase. However, since we do not 
include thermal fluctuation (noise) in our kinetic equation (|2]l, 
the a -/3 front is virtually arrested as there exists a small barrier 
between a well and /3 well (Fig.[Tll. 




Figure 7: A composite nucleus of stable y stripe (white) sandwiched by a less 
stable a phase (black) embedded in most metastable /3 phase (gray) at / = 200. 
This configurations is kinetically arrested and cannot grow even though there is 
a seed of the most stable y phase (white). 



The appearance of this kinetically arrested metastable /3 
phase around the composite nucleus made from the most sta- 
ble y core wrapped by the least stable a skirt is not due to the 
planer symmetry. Figure [8] shows the evolution of the multiple 
of circular composite nuclei consist of most stable y phase of 
radius 5 surrounded by a thin layer of least stable a phase of 
radius 8 embedded in a least stable /3 phase. Again, these com- 
posite nuclei are kinetically arrested and cannot grow. Only a 
fusion of two colliding nuclei occurs near the center. The re- 
sulting oval nucleus of y wrapped by a thin layer of metastable 
a is static and stable again. Then the metastable /3 phase envi- 
ronment survives even though there are seeds of most stable y 
phase. This configuration called boiled-egg structure was pre- 
dicted from experiments , 27 , 28l] and was found by the cell- 
dynamics simulation of three-phase system by the author |01 
for the special case of single composite nucleus and the a-y 
equilibrium. 




Figure 8: Composite nuclei with the most stable circular core of y phase sur- 
rounded by a thin layer of less stable a phase embedded in least stable /? phase. 
Only a fusion of two colliding nucleus is observed near the center. This boiled- 
egg configuration is stable and kinetically arrested. 



Since we start from the static configuration, not only the evo- 
lution of the most stable y phase by a — > y transition is prohib- 
ited but also the evolution of the surrounding a phase hy /3 a 
transition cannot be initiated spontaneously as the free energy 
barrier between a phase and /3 phase exists (Fig.[TJ. Naturally, 
the introduction of the thermal noise initiates /3 —> a transi- 
tion 101 . However, this kinetically arrested phase implies fairly 
long lifetime of these metastable configuration 



that is characterized by a single non-conserved order parame- 
ter. We have considered three cases when a relative stability of 
three phases, a, /3 and y phases are diff'erent. We have found 
several interesting scenarios of the kinetic of phase transition in 
a three-phase system: (a) the appearance of stable third phase 
from the interface of two metastable phase, (b) the appearance 
of dynamically stable third metastable phase from the moving 
interface between stable and metastable phase, and (c) the ap- 
pearance of kinetically arrested third metastable phase around 
the composite nuclei made from the stable phase surrounded by 
a thin layer of metastable phase. Although, the last two exam- 
ples have already been predicted theoretically fT?, TT], or nu- 
merically liTJ for special configuration, we showed in this paper 
using cell-dynamics simulation that these predictions are qual- 
itative correct even for multiple nuclei in two-dimensions and 
for more wide variety of free energy landscape. 
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4. Conclusion 

In this paper, we have used the cell dynamics method to ex- 
plore various scenario of phase transition in three-phase system 
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